Symmetry restoration at finite temperature with weak magnetic fields 
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We study symmetry restoration at finite temperature in the standard model during tlie electroweak 
phase transition in the presence of a weak magnetic field. We compute the finite temperature 
effective potential up to the contribution of ring diagrams, using the broken phase degrees of freedom, 
and keep track of the gauge parameter dependence of the results. We show that under these 
conditions, the phase transition becomes stronger first order. 
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I. INTRODUCTION 

The properties of the universe's primordial plasma at 
high temperature play an important role in the attempts 
to explain several outstanding questions in eosmology 
such as the origin of the matter-antimatter asymmetry. 
Starting from symmetric conditions, it was shown by 
Sakharov [l[ that three ingredients are needed to develop 
a baryon asymmetry: (i) Baryon number violation pro- 
cesses, {a) C and CP violation and {iii) departure from 
thermal equilibrium. These conditions are met in the 
Minimal Standard Model (MSM) during the Electroweak 
Phase Transition (EWPT) provided that this be first or- 
der. It is however well known that neither the amount 
of CP violation nor the strength of the phase transition 
are enough to produce and subsequently preserve a pos- 
sible baryon number at this stage of the universe evolu- 
tion [2„i3|]. 

A great deal of effort has been devoted to put forward 
viable scenarios that increase the amount of CP violation 
and/or make the EWPT stronger first order 0, Q . At 
the same time, in the last years, the study of the effects 
of primordial magnetic fields on other cosmological pro- 
cesses @-l3 , have also been given attention, including the 
very same EWPT %^. 

Since magnetic fields before the EWPT belong to the 
t/y(l) group, this kind of fields properly receive the name 
of hypermagnetic. The expectation is that in their pres- 
ence, the order of the EWPT is increased, in analogy with 
the case of superconductivity, where an external mag- 
netic field changes the order of the phase transition from 
second to first due to the Meissner effect [3, [TH [12] ■ 

The description of the physical processes able to gen- 
erate magnetic fields is an old problem in cosmology [Ij] . 
The general approach is to identify mechanisms for the 
generation of seed fields [ij, [l^ which can later be am- 
plified into fields on larger scales. For instance, certain 



types of inflationary models can pro duce magnetic fields 
extending over horizon distances J16l |. Other possibilities 
make use of cosmological phase transitions [13 through 
formation of dipole charge layers on surfaces of phase 
transition bubble walls. Although at present there is no 
conclusive evidence about the origin of magnetic fields, 
their existence prior to the EWPT cannot certainly be 
ruled out. They have been observed in galaxies, clusters, 
intracluster medium and high redshift objects [18| . 

The presence of magnetic fields in the early universe 
can leave its imprint in a variety of phenomena. Mag- 
netic fields directly interact with ionized baryons and 
indirectly influence photons, through the tight coupling 
between baryons and photons at that epoch. Cold dark 
matter is also indirectly affected, through gravitational 
interaction. Thus, a primordial magnetic field imprint 
can be searched for (through temperature anisotropics, 
statistics and polarization) in the Cosmic Microwave 
Background (CMB), formation of large scale structure, 
gravitational wave background and the nucleosynthesis 
processes. A homogeneous magnetic field would give 
rise to a dipole anisotropy in the background radiation. 
On this basis. Cosmic Background Explorer (COBE) re- 
sults give an upper bound on the present equivalent field 
strength of Bq < IQ-^ G [11 [ll|. A lot of work has been 
done recently in the field of the observational constraints, 
considering different aspects of the interaction of primor- 
dial magnetic fields with the CMB, as well as different 
features and scales of these cosmic fields (see for exam- 
ple 12 0| and references therein). In particular, in a recent 
work [2l[ , using both CMB and large-scale structure ob- 
servational data, the authors place new constraints on 
primordial magnetic fields, both on their magnitude and 
power spectral index. This last parameter is important 
for the discrimination between models of magnetogene- 
sis. For all cosmological parameters, the same priors as 
those adopted in the WMAP analysis arc used. Standard 



cosiTLological parameters are mainly constrained by low 
multipoles, while primordial magnetic fields dominate on 
small angular scales (high multipoles). The degeneracy 
among magnetic fields parameters is broken by the dif- 
ferent effects of primordial magnetic fields on the matter 
power and CMB radiation spectra. An interesting result 
of this analysis is that the probability distribution of the 
parameters has a maximum for nonzero values of both 
the magnitude and the power spectral index. Although 
these values are still consistent with zero magnetic fields 
(thus implying only upper limits on the field strength of 
~ 10~^ Gauss on a comoving scale of 1 Mpc) they suggest 
the possibility of a detection with forthcoming data. 

Recall that the development of the EWPT can be de- 
scribed by means of the MSM effective potential (EP) 
where the order parameter is the Higgs vacuum expecta- 
tion value V [22, [23] . Starting from a very high tempera- 
ture, the minimum of the theory happens for v = 0, the 
so called symmetric phase. Decreasing the temperature 
produces the EP to develop a secondary minimum, in 
case the phase transition is first order. Since the phase 
transition does not start until the secondary minimum 
is degenerate with the original one, its development in 
the high temperature phase can be described in terms 
of MSM symmetry restored degrees of freedom. Parti- 
cle masses depend parametrically on the order parame- 
ter and acquire their fixed values only when the phase 
transition is completed. 

Within the particle physics calculational methods, it 
is common to work using the degrees of freedom in the 
symmetry restored phase. A drawback of this scheme is 
that, since the gauge boson mass matrix is not diagonal, 
contributions from gauge bosons to the EP are harder 
to compute. Another general aspect of calculations in- 
volving MSM degrees of freedom is that they are usually 
performed for a given value of the (covariant) gauge pa- 
rameter, which makes it impossible to keep track of a 
gauge parameter presence, if any, in phase transition re- 
lated obscrvables. 

In a recent work [2J] , the development of the EWPT 
in the MSM from the symmetric phase has been studied 
in the presence of a weak external hypcrmagnetic field 
up to the contribution of ring diagrams. The main result 
from that work is that the presence of the field strength- 
ens the first order nature of the phase transition. The 
result is in agreement with calculations performed at a 
classical [il[ and one- loop levels [i2[, as well as with lat- 
tice simulations Q. On the other hand, other analytical 
non-perturbative approaches where the MSM finite tem- 
perature EP is studied for the case of strong magnetic 
fields [2^, [2g , reach the conclusion that these fields in- 
hibit the first order phase transition. These works at- 
tribute the result to the contribution of light fermion 
masses which are generally neglected in other computa- 
tions. However, they also neglect to consider the infrared 
cutoff provided by thermal masses, which casts doubts on 
their conclusions. It is thus important to study the phase 
transition using the broken phase degrees of freedom. 



since it will help to set up the stage for a detailed anal- 
ysis of the phase equilibrium conditions as the EWPT 
develops. 

In this work wc generalize the analysis of the develop- 
ment of the EWPT in the presence of a constant magnetic 
field, working with the degrees of freedom in the broken 
symmetry phase where the symmetry of the theory has 
been reduced from SU{2)l x U{1)y to t/(l)em- We work 
with an arbitrary value of the (covariant) gauge param- 
eter in the weak field limit and up to the contribution of 
the ring diagrams, that have been shown to be crucial 
for the description of the long wavelength properties of 
the theory ^^. We find that there is a small gauge pa- 
rameter dependence on observables such as the critical 
temperature and the position of the broken phase mini- 
mum, but otherwise confirm that the presence of a weak 
magnetic field, increases the order of the phase transition. 

Also, as we will see in the next section, we work ex- 
plicitly with the assumption that the hierarchy of scales 
eB <^ m? -^ r^ is obeyed, where we consider m as a 
generic mass of the problem at the electroweak scale. 
Since at the phase transition m '^ v, the effective cou- 
pling in the perturbative expansion becomes g^T/v and 
thus when v/T is small {T/v is large), the perturbative 
series is not reliable. This issue has been discussed in 
detail by Arnold and Espinosa [2^ where they use a 
power counting argument whereby the leading tempera- 
ture terms (after accounting for the thermal contribution 
to the effective mass) give rise to corrections proportional 
to g^T/m. Notice however that in this work we do not 
attempt to claim that the hierarchy of scales that we 
worked with makes v/T grow to the extent that both, the 
sphaleron erasure bound is avoided nor the perturbative 
expansion problem is solved. Rather, our purpose is to 
explore the effective potential method to see if a hint of 
a growth of v/T is possible. As we will see, such growth 
does happen and in that sense the results are encourag- 
ing. Another recent study of the effective Lagrangian in 
the MSM using broken phase degrees of freedom at finite 
temperature and density can be found in Ref. [291 . 

The paper is organized as follows: In Sec. |TT]we write 
down the MSM using the degrees of freedom in the bro- 
ken symmetry phase. In Sec. lIIII we lay down the formal- 
ism to include weak magnetic fields in the computation 
of charged particle propagators. In Sec. lIVi we work with 
these degrees of freedom to compute particle self-energies 
that are used in Sec. |V] to compute the MSM EP up to 
the contributions of ring diagrams. In Sec. IVII we study 
this EP as a function of the Higgs vacuum expectation 
value and show that the order of the EWPT becomes 
stronger first order in the presence of the magnetic field. 
We also study the gauge parameter dependence of the 
phase transition parameters. Finally, we conclude and 
discuss our results in Sec. IVIII and leave for the appen- 
dices the listing of the one-loop MSM self-energies in the 
presence of weak magnetic fields, for arbitrary values of 
the gauge parameter. 



II. MINIMAL STANDARD MODEL 

Quantization of the MSM in an arbitrary covariant 
gauge using the Faddeev-Popov technique involves not 
only physical fields such as the gauge bosons W^, Z^, A^, 
the neutral Higgs h, and fermions, but also charged and 
neutral Goldstone bosons x^iXSi respectively, which are 
fields that do not propagate asymptotically. Neverthe- 
less, these last fields appear in loops. Working in a 
renormalizable gauge of the t'Hooft type, the so called 
i?j gauge, in order to compensate contributions from the 
spurious components of the gauge fields, it is necessary 
to also introduce Faddeev-Popov ghosts for each of the 
gauge boson fields rj^^^-qz and 77^. Thus, the Lagrangian 
also depends on the choice of the gauge fixing term. 

In order to consider all the contributions to the MSM 
EP, we write the complete Lagrangian sector by sector, 
in an arbitrary R^ gauge, after the symmetry breaking. 



C — Ch + Cgb + Cf + Cy + Cgf + Cfp, 



(1) 



where Ch, Cgb, Cf, Cy, Cgf and Cpp are the La- 
grangians for the Higgs, gauge boson, fermion, Yukawa, 
gauge fixing and Faddeev-Popov sectors, respectively. 
The Lagrangian for the Higgs field is 



Ch = {Df,'^y{D^'<S>) + c^($T$) _ A($T$) 



(2) 



The SU{2)l x U{1)y symmetry of the Lagrangian in 
Eq. ([T]) is broken spontaneously by introducing an ad- 
equate Higgs scalar $, which develops a vacuum expec- 
tation value V. The most economical choice for 4> is an 
SU{2)l doublet of complex fields with hypercharge (y) 
+1, given by 



V2 



X2 + iXl 



(3) 



We take h as the physical Higgs field with a mass given 

by 



'ff 



3Aw" 



(4) 



From the kinetic energy term on Eq. ^ one obtains 
the gauge boson masses by finding the mass eigenstates 
in the broken symmetry phase, which leads to 



i-w 



gV/4, m| = (g2+.g'2)z;V4 (5) 

and the relations between the gauge fields W^,B^ and 
W^ , Z^j, A^ before and after symmetry breaking, respec- 
tively 



w^ = ^{K^'K) 



z„. 



V? 



A,^ 



V? 



{gWl - g'B^) 
{g'Wl+gB^) 



(6) 



Hereafter g and g' are the SU{2)l and U{1)y gauge cou- 
plings, respectively. 

In order to obtain the interaction terms between the 
Higgs field and the gauge bosons, it is convenient to 
rewrite the covariant derivative 



Dn 



^^ + W^W^+^^B' 



(7) 



in terms of the mass eigenstates where B' = B^ 



BT 



This is explicitly accomplished by applying D^ in terms 
of the states defined in Eq. ([6]), to the Higgs field in 
Eq. ([3]). On the other hand, the self-interactions of the 
Higgs field, are obtained from the scalar potential in 
Eq. ©. 

The kinetic energy from the SU{2)l and U{1)y gauge 
bosons is 



I 



f 



where 






dnW, - a„W„ - .gW„ X W„ 



(8) 



B'nv - 9^Bl - di,B' 



(9) 



After symmetry breaking and introduction of the mass 
eigenstates, one obtains from Eq. ([8]) the kinetic term as 
well as interaction terms involving three and four gauge 
fields. 

The Lagrangian for the fermion sector is 



Cf =*fl^*fi + fi^^L, 



(10) 



where ^l,r = ^(iTTs)*- Once again, the fermion-gauge 
boson interactions are obtained by writing the covariant 
derivative in terms of the mass eigenstates after symme- 
try breaking. 

We work in the limit where all fermion masses are neg- 
ligible, except the top quark mass. Therefore, the main 
contribution to the Yukawa sector is 



Cy = VtqL^tji + h.c, 



(11) 



where yt is the top quark Yukawa coupling, q^ is the third 
family quark doublet, tp is the right handed top quark, 
$ = ia2^* and CT2 is the second Pauli matrix. After sym- 
metry breaking this Lagrangian yields for the top quark 
To determine the interaction terms 



mass mt 



Vt 



between fermions and the Higgs field, it is necessary to 
rotate the fermions to the mass eigenbasis. 

In the R^ gauge, the gauge fixing Lagrangian is 



C 



af 



1 

2? 
1 



I 



{d^Wl - -igvx') 



i\2 



1 



^^{d''Z,-:^£,g'vx,f 



(12) 



where i — 1,2 and ^ is the gauge parameter (which for 
simphcity we set it to be the same for all gauge bosons) 
and rrii stands for the Goldstone boson masses, given by 



Xv' 



c' + ^g' 



and 



\ 2 2 

\v — C 



/2n« 



+a9'+9') 



4' 



(13) 



(14) 



respectively. 

Since we work with an arbitrary value of the gauge 
fixing parameter ^, the ghost fields acquire masses pro- 
portional to the corresponding gauge boson ones nigb 



"rigb 



^771 



gb- 



(15) 



Notice that these fields do not propagate asymptotically, 
thus appear only in internal Feynman diagram lines. 

In order to find the corresponding Faddeev-Popov 
ghost interactions, the ghost fields are rotated like the 
gauge fields in Eq. ^. We emphasize that by working 
with an arbitrary value of ^, the ghost fields do con- 
tribute to the w-dependent part of the one-loop EP. Note 
that the EP is in principle a gauge dependent object j30( . 
however, physical quantities obtained from it should be 
gauge independent j31| . 



III. CHARGED PARTICLE PROPAGATORS IN 
THE PRESENCE OF A MAGNETIC FIELD 

We work with MSM degrees of freedom in the sym- 
metry broken phase, where the external magnetic field 
belongs to the C/(l)em group. To include the effect 
of the external field, we use Schwinger's proper time 
method [32[. In the broken phase, we have three kinds 
of charged particles that couple to the external field, 
namely: scalars, fermions and gauge bosons, whose prop- 
agators are 



Db{x,x') = (f>{x,x') 



Sb{x,x') = 4>ix,x') 



G^{x,x') — (j){x,x') 



(2^ 



-ik-{x—x ) 



DB{k), (16) 



e-*(^-^')S'B(fc), (17) 



-*(^-^')G^''(fc),(18) 



respectively. The phase factor 0(x, x'), that breaks trans- 
lation invariance, is given by 



'^ = p«=L"' d«^[sr+I^M.(«-^')1 



(x,x ) = e 



(19) 



where the vector potential 5°"' = -j(0, y, — x, 0) gives 
rise to a constant magnetic field of strength B along the 
z axis and Fj^J ~ 9^5°"' — d^B°^^ is the external field 
strength tensor. 



The momentum dependent functions Dsik), Ssik) 
and G^ (k) are given by 



iDeik) 



ds 



cos eBs 



tan eBs 



X exp ■{ is(kf: — k'\ ttt-^ + ie) 

' " eBs 



(20) 



iSsik) 



ds 



cos eBs 



X exp ^ is{kn — k^ 



taneBs 



eBs 



{nif + |(f||)e 



ieBsa^ 



^A 



and 



^G^^ik) = 



oo 



ds 



cos eBs 



2 _t,2 taneBa 



(21) 






cos eBs 



'9':r+ie"''") 



2eFs\p.i' 



,-is( 



™36-»eJ 



= -*s(C' 



m„.-ie] 



'gb 



X [{k^' + kxF^'^{{i&u{eBs))/B)) 
X (k" + kpFP''((t&n(eBs))/B)) 

_ it (F^'^ + g>'^Bitin{eBs))\ \ . 



(22) 



We use the metric tensor g^'^ = diag(l, — 1, — 1, — 1) in 
which g^^^ = g';^." — g^ and the notation fcn ~ k^ ~ k^, 

Y^jij^. We also use 



fci 



kf + k2, and cr'^ = 17^7^ 



(g2ef .)M. = g^^^ cosi2eBs) - F^- El^lifM (23) 

B 

where throughout u^ and 6^ are four-vectors describing 
the plasma rest frame and the direction of the magnetic 
field, respectively. In the rest frame, these are given by 



" = (1,0,0,0), 



(0,0,0,1). 



(24) 



It has been shown that, by deforming the contour of 
integration, Eqs. pop and ((2T|) can be written as [3^, |3J| 



(25) 



.!>.(.) = 2. £ (-l)'i'(^)-^ 



k?, - (2/ + l)eB-m2 + ie' 



/=o '"11 

^ ' ^ k^\- 2ZeB - mj + le k]_ 
where di{a) = (-l)"e-"L-i(2a), <, = ddn/da, 

"i^-fc|l 



(26) 



D = {nif + }i\\) + }i^- 



kl 



(27) 



and Li, ij" arc Laguerrc and Associated Lagucrrc poly- 
nomials, respectively. It can also be shown that by per- 
forming a similar analysis as in Refs. [33, [SJ], Eq. (|22] 
can be written as ISSlI 



1 oo 



^G'Bik)= EEt^ 



i„-^ 



2i(-l)'e-^ 



A=-l i=0 II 



kf. - {21 + 2X + l)eB - ml^ + ie 



rptlV 



"^'b(i - OP""" 



fc| - {21 + 2A + l)eB - em2^ + ze 



-§ 



where 



r^- = E [.9'^"^ (1^1 - 1) 



A=-l 






Ff"" / X 



m' 



ptiv 



J_[(fcMfc.„,J^A<.) + ('^-A 



B \4i. 



m 



gb 



B 



kpkf^FP-^ eB d 



+ ( 



B 

kxkpFf'^FP'' 



B^ 



){teBf 



d{klr' 



(28) 



(29) 



In order to determine the appropriate order of energy 
scales during the development of the the EWPT, one re- 
sorts to bounds on the strength of the magnetic fields 
imposed by cosmological processes in the early universe. 

I 



The relation between the strength of large scale mag- 
netic fields and temperature is obtained from the re- 
quirement that the magnetic energy density pmag ~ B^ 
should be smaller than the overall radiation energy den- 
sity Prad ^ T^ at nucleosynthesis, in order to preserve the 
estimated abundances of light elements. With this, one 
obtains the simple bound B < T^ [3^|. Furthermore, to 
guarantee stability conditions against the formation of a 
VF-condensate [33 , one obtains that the field is also weak 
compared to rriyy. We want to study the EWPT transi- 
tion via the analysis of the EP when its minima change 
in the presence of magnetic fields, we will work explicitly 
with the assumption that the hierarchy of scales 



eB <^m'^ -^ T^ 



(30) 



is obeyed, where we consider m as a generic mass of the 
problem at the electroweak scale. It is important to ad- 
dress the same problem with the hierarchy of the mass 
and the magnetic field scales switched and this will be 
done elsewhere [3^]. We can thus perform a weak field 
expansion in Eqs. ([25]) . ([26]) and ([28]) which allows us to 
carry out the summation over Landau levels to write the 
scalar, fermion and gauge boson propagators as power 
series in eB, that up to order (eS)^ read as [33l - [35l | 



DB{k) 



1 



fc2- 



1 



{eB)' 



2{eBykl 



(31) 



^"'^^r-m? 



^ + nif 75^^(/c|| + mf){eB) 



(fc2 



TO?)-^ 



2{eBfkl 



(fc2 



m]Y 



{mj + /(f|| + /(fj 



A-2 



kl 



and 



), (32) 



G%^{k) = ^ 




i{eB) 



2 f kxF^^knFP- 



B' 



(fc^-"^s&) (^^-^"^ffb) y 



(33) 



respectively. Notice that for the assumed hierarchy of en- 
ergy scales to be valid, one should stay away from work- 
ing with the value for the gauge parameter ^ = 0, since, 
according to Eq. ([15]) this choice would lead to the van- 
ishing of the ghosts masses. 



IV. SELF-ENERGIES 

In this section we compute the MSM sclf-encrgics that 
are in turn used for the computation of the ring diagrams 
in the EP. 

It is well known that in the absence of an external 
magnetic field, the MSM thermal self-energies are gauge 
independent when considering only the leading contri- 
butions in temperature |39|. However, as we will show, 
when considering the effects of a weak external magnetic 
field, these sclf-encrgics turn out to be gauge dependent. 

In what follows, we work in the imaginary-time for- 
malism of thermal field theory. First, we note that the 
integration over four-momenta is carried out in Euclidean 
space with fcg = ik4, this means that 

7^-^/^- (34) 

(27r)4 J (27r)4 ^ ' 

Next, we recall that boson energies take on discrete val- 
ues, namely k^i — a;„ — 2mrT with n an integer, and 
thus 



d^kE 

(2^ 



^E 



_fk_ 

(2^ 



(35) 



A. Higgs boson 

Figure [1] shows the diagrams that contribute to the 
Higgs boson self-energies affected by the magnetic field. 
Let us explicitly compute the momentum independent di- 
agram shown in Fig. [ija) for a single scalar field running 
in the loop, H^^'^. In the weak field limit, its expression 
is 

r^£'B(w„,k;m^), (36) 



H 



H-S 



A V-- 
4^E 



(27r)3 



where Db is given by Eq. ([3T|) . 

Using the Euclidean version of Eq. ((3T|). we have 

4"i (2^) 



H 



H-S 



1- 



(2^)3(c.2+k2 



m2) 



(w2+k2+m2) 



2^2 



2kl{eB)^ 



(37) 



(w2 +k2 +?7l2)3^ 

The integrand in Eq. p7p contains terms whose general 
form is 

/.(k;m^)^ , . J^ „, . (38) 




S : 


Higgs and Goldstone bosons 


G 


Z and W gauge bosons 


F 


fermions 


Gh 


: ghosts 



FIG. 1. Feynman diagrams for the Higgs boson self-energy. 



The integral involving the sum over Matsubara frequen- 
cies is performed by resorting to Refs. [40, |41|. For the 



terms with n ^ the result is 



^E/(0zk-.-/.(k;..^) 



(2r) 



(-2_^^^d+2a+2t-2Q "*- V 2 ' "^z ,,2e 



r(i) 



(47r)''/2r(a) 



m \2i 



whereas for the term n = we get 



(39) 



d'^k k2 



T r(| + a) 



(27r)<' (k2 + m2)" (47r)''/2 r(f) 
r(a-|-a) / 1 ^"---« 



r(a) \m 



(40) 



where C is the modified Riemann Zeta function, fi is the 
energy scale of dimensional regularization and d ~ 3 — 2e. 

In terms of Eq. ^E^, Eq. ^7^ is written as 



/i(k;m2) 



n 

(e,B)2 



^2^ 



/3(k;mf)-2k^/4(k;TOf) . (41) 



Using Eq. (p9|) for the terms with n 7^ and Eq. (|40)) for 
n = 0, and after considering the contributions from all 
scalar fields running in the loop, we get 



rp2 



1 



2ttT 



(2toi + 1713, + 37714) 



-A 






T , C(3)m2 



?7ii Att^T^ 



(42) 



where tti^ stands for the Goldstone boson masses, given 
by Eq. Ipj). 

The contribution from the diagram in Fig. W[b) for a 
W boson running in the loop is 



tH-w _ 9 
^'' ~ 2 



nf-"^ =^T^ 



d^k 



(27r)3^ 



C^ {ujn,'k;mw)- (43) 



We proceed to compute Eq. ^^ , as well as the rest of 
the self-energies, in the infrared limit [ij (70 = 0,q ^- 0, 
unless explicitly stated. 

In terms of the function /„ defined in Eq. psp we can 



write Eq. (gSj) as 

„2 



n 



H-W 



2™^ 



T 



4^7 (2^)3 I 



4m^/i(k;m 



w) 






6(/2(k; m^) - /2(k; ^m^)) 

- {3loI + 7k2)(/3(k;m2^) - hiK^ml^)) 
+ ik^cjl + k2)(/4(k; ml.) - /4(k; ^m^^)) 

+ 36TO^/3(k;TO^) + 16k^m^/4(k;TO^) |. 

(44) 

Using Eq. dSHl) for the terms with n 7^ and Eq. (HDl) 
for the terms with n = and including all contributions 
from the gauge bosons in the loop, we get 
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ttT 

(eS)V 
647rm2^, 



{2g^mw + (g^ + g'^)mz) 

35 \ T , llC(3)m2,. 



3^1/2 



mw 



37r3r2 



(45) 



In a similar fashion the contributions from the diagrams 
(c) to {g) in Fig.[Tl in the infrared limit, can be computed 
and are given Appendix A for an arbitrary value of the 
gauge parameter ^. 

The leading contribution to the Higgs self-energy is 
obtained by adding expressions (a)-(g) in Appendix A, 
keeping the leading term for each contribution. The 
result is gauge parameter independent and is explicitly 
given by 



Hi 



4 



rff 



2X + f 



(46) 



Notice that for large values of the top quark mass, namely 
a large coupling constant /, a perturbative calculation 
is not entirely justified. Nevertheless, here we consider 
our calculation as an analytical tool to explore this non- 
perturbative domain. 

In expressions (a)-(g) for the Higgs self-energy in Ap- 
pendix A, we have kept terms representing the lead- 
ing contribution of each kind arising in the calculation, 
namely, terms that after extracting a T^ factor, are of or- 
der (eB)VT*, m/T, v^/Tm and {eB)^/Tm^, where m is 
a generic mass. For the hierarchy of scales considered, the 
first kind of terms can be safely neglected. Also, terms of 
order m/T are small. These kind of terms include ratios 
of gauge boson masses to the temperature, which accord- 
ing to Eq. (O leads to ratios of v/T. When the masses 
appear in the denominator and are the scalar ones, the 
terms are potentially dangerous since their square can 
become negative. However, as we will soon show, these 



kind of terms arc naturally canceled in the effective po- 
tential and substituted by terms where the mass that 
contributes is the thermal one, defined by 



n,. 



(47) 



When the gauge boson masses appear in the denomina- 
tors throughout, one should recall that such masses are 
to be thought of as being computed in the broken phase. 
However, since our analysis considers the effective poten- 
tial as a function oi v > 0, we will implement such re- 
striction by replacing in the effective potential the gauge 
boson masses appearing in denominators by the thermal 
ones, whose definition, in analogy to the scalar case, is 
given as the square root of the sum of the gauge boson 
mass squared and the leading term for the corresponding 
self-energy, to be defined shortly [sec Eqs. (p^ and (p5|) ]. 




[c) 



\ 



y S 



B. Gauge bosons 

To express the gauge boson self-energies, in the pres- 
ence of the external field, we notice that we have three 
independent vectors to our disposal to form tensor struc- 
tures transverse to the gauge boson momentum g^. 



namely q'' 



and &'', where these last two vectors 



are given in the rest frame by Eq. ([M]) . This means 
that in general, these self-energies can be written as lin- 
ear combinations of nine independent structures |43l . |44| . 
Since we are interested in considering the infrared limit, 
qo = 0, q — J- 0, only u^ and b^ remain. Notice that 
the correct symmetry property for the self-energy is 
n^''(q) = IT^''(— g) [ill. However, in the infrared limit, 
this condition means that the self-energy must be sym- 
metric under the exchange of the Lorentz indices and 
therefore we can write: 

H^" = H^Q^"^ + U"R^"' + H^^^'' + H^^g'^", (48) 



where 



IM!^ 



/A',,'^ 



S^'"' :^u^'h'' + u''b'', (49) 

and the transversality condition g^H'^" = is trivially 
satisfied in the infrared limit. Also, note that working in 
the rest frame of the medium [see Eq. ((24)) ] , 



noo = n« + H 



M 



(50) 



Figure [5] shows the gauge boson self-energy diagrams, 
which involve Goldstone bosons, fermions as well as 
gauge bosons in the loop. Let us explicitly compute 
the component H^'' for the self-energy diagram shown 
in Fig. [5|a) for the case where the external particle and 
the loop particle are VK-bosons. The explicit expression 
for all the components of the tensor is 

^^^- >^(g) = ig\2g^P9''- - g'^^g"- - g'^'^ g^'^) 

^GB,Ak)- (51) 



H 



(a) 




S : Higgs and Goldstone bosons 
G : 7, Z and W gauge bosons 
F : fermions 
Gh : ghosts 



(27r) 



FIG. 2. Feynman diagrams for the Gauge bosons self-energies. 



Using the Euclidean version of the gauge boson propa- 
gator G'^ obtained from Eq. (|33l) . and in terms of the 



function /„ defined in Eq. ([38|) , we can write the compo- 
nent that is being considered at finite temperature as 



n. 



w-w 
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3m^ 



^E 






|9m^-/i(k;TO^) 



3k2(/i(k;m2,.)-/i(k;Cm2^)) 

(eB)2 [-3(/2(k;774)-/2(k,0;^m^) 

7k2(/3(k;m2^,)-/3(k;Cm2^,) 

4k'*f/4(k;m2^)-/4(k;Cm2^) 



15m^/3(k;m^) + 12k^m^/4(k; to^)] |, 



(52) 



where hereafter 11''° ^^ = II^'' for gauge boson self- 
energies. 

Using Eqs. ([39]) and gOl) into Eq. ^, including aU 
contributions from gauge bosons in the loop and keeping 
only the leading terms we get 



n! 



2g2r2(^ + l)- 



(2 + e^/2)g2r 
47r(g2+g'2) 



g2(ei?)2 ^ 4r f^ 2S_\ ^ 7a^)m'w 

^1/2 



USnm'^- 






7^37^2 



(53) 



We now turn to compute the expressions for 
n^ depicted in Fig. [2l^b). For the case where 

the external particle is a W^-boson and the internal ones 
are a neutral gauge boson and a charged Goldstone bo- 
son, its explicit expression is 



nr 



(9) 



ah''' 



n 



+ 9' 



[m 



2 
w 



mt 



(27r) 



jg^,pG'"'ik)g,^DB{k~q), (54) 



where the upper (lower) case corresponds to an internal 
Z (7)-line. Notice that since the net charge fiowing in the 
loop is not zero, the phase factor referred to in Eq. (fT9|) 
does not vanish and should in principle be also included 
in the computation of this self-energy. Nevertheless, as 
we are interested in this expression to compute the con- 
tribution of the ring diagrams to the EP, that in turn is 
represented by closed loops, the phase factor becomes the 
identity and therefore it does not need to be computed 
for individual self-energies. The same applies for the rest 
of the self-energy diagrams where there is a net charge 
flowing in the loop and consequently, we do not consider 
computing their corresponding phases. 

Using the Euclidean version of Eq. ([5T]) and Eq. ([551) , 



we have 



n, 



W-Z('y)-S_ 



g'g" 



(m| - m^) 



g2 + g,2 

^E 

n 

|m|(^)/ii(k,0;m,^m| 
2 



'w 



-w: i^lii{k,0]mf,nt.z(^j) 



_fk_ 

(27r)3 



(7)^ 



^ mi,„,^) -/ii(k,0;mj^^m|(^) 



|_4k TO2(^)/4i(k,0;m,,m2(^)) 
-3TO|(^)/3i(k, 0; m,^ to|(^)) 

-3w,^j(/3i(k, 0; m- , m|(^)) - /3i(k, 0; m- , ^m|(^-|)) 
+4a;2 k^ (/41 (k, 0; m^ , m|(^) ) - /41 (k, 0;ml im^i) ) ) } • 

(55) 



The integrand in Eq. (|55p contains terms whose general 
form is 



/Q/^(k,q;m , ^m ) 



1 



[u;2 +k2 + m2]" 

1 
[a;2 + (k-q)"+em2]/3' 



(56) 



We make use of the Feynman parametrization to write 
laji as 



/a^(k,q;TO ,^m ) 



2,2. r(« + /3) 



r(a)r(/?) 

[c^2+k'2(a;)+m'2(x)]"+'3' 



(57) 



where 



k'(x) = k - (1 - x)q 
m'^{x) = 77i^(^ + x - ix) + a;(l - a;)q^ (58) 



and r is the Gamma function. Notice that according 
to our findings in Ref. 41|, this parametrization is al- 
lowed since we are pursuing a description in the infrared 
limit. To carry out the sum over Matsubara frequencies 
together with the integration in Eq. (|55p. onc e ag ain we 
resort to the results in Ref. [40] (see also Ref. [2J]). The 
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explicit result, generalized to include the fermion case is 

(27rT)''+2°+2*-2("+/5) r(| + a) 



T{a)T{(5) 



Til) "^ 



xr(j + a + /3---a) 

X l"dx.-'il-xf-'('4^y\ (59) 







27rT 



where for fermions the sum runs over all integers and 
Z = 1/2. For the terms involving the n = Matsubara 
frequency for bosons, we use the result 



r(| + a)r(a + /j-|) 

r(f) r(«)r(/3) 



T 



(47r)d/2 



„1 / 1 \ 2a+2/3-(i+2a 

X / dx x""i (1 - xf-^ ( — — 
^ \m'{x) 



(60) 



Using Eq. ([591) for the terms with n ^ and Eq. (jeOj) 
for n = 0, and after considering the contributions from 
all gauge boson fields running in the loop, we get 
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77ii m\ 
m,^ 



32 ((72 + (7'2) 7rr77,2(TOx + TO^)3 y TOl 



(ei3)277^2^C(5) 
20487r6r4 



y^- 155^ -145V 
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(61) 



The rest of the diagrams depicted in Fig. [2] can be com- 
puted in the same manner and the result is given in Ap- 
pendix B. For the expressions corresponding to diagrams 
(a) - [g) in Fig.[2J some remarks are in order: First, no- 
tice that terms of order rUgh/T, where mgb stands for any 
of the gauge boson masses, are proportional to v/T and 
thus, as in the discussion of the scalar self-energy, these 
terms are small when the position of the minimum in the 
broken phase is small compared to the critical tempera- 
ture. Second, when the gauge boson masses appear in the 
denominators throughout, we should keep in mind that 
the analysis is valid in the broken phase, where the vac- 
uum expectation value of the Higgs field is different from 
zero. However, since the analysis considers the effective 
potential as a function of w > 0, we will implement this 
restriction by substituting the gauge boson mass by the 
thermal ones, when this mass appears in denominators. 
This is meant to regulate the singularity at u = and 
does not make a numerical difference while allowing us to 
treat the effective potential as a function in the domain 
V > 0. The gauge boson thermal masses are defined by 



mw 



mz 



'"W 



(n' 



VKA: 



(n^) 



zJl; 



.^-y'(nV)i, (62) 

where (n^^)i, (H^)! and (n^)i are the leading terms 
of the gauge boson self-energies made up by adding the 
corresponding terms in Appendix B for each gauge bo- 
son. These leading terms turn out to be gauge parameter 
independent and are explicitly given by 

_ 11 



(n^)i = ^-^9'T^ 



Q. _ll(/+5'^) 



(uQ] - ^^ g'g'' t2 



(63) 



As is described in Ref. 133|, the other non zero compo- 
nents of the gauge boson self-energy are negligible 



n" = n22, 033^0(7772) 
n03 = n^o ~ oieB). 



(64) 

From Eqs. (gS]) and ^U^, we see that since 11*'^ = -II" 
thenn™ ~nQ, and thus 



U^"' ~ n^Q'"''. 



EFFECTIVE POTENTIAL 



(65) 



A. One-loop 



In the standard model the tree level potential is 



(66) 
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To one loop, the EP receives contributions from each 
sector, namely 

where in general each one of these contributions is given 

by 



v^^Hv) 



T 
2 



j:J^^^MDi-n,k)-^]),m 



with D standing for either the scalar, fermion, gauge bo- 
son or ghost propagator, and the trace is taken over all 
internal indices. We first discuss the thermal contribu- 
tions coming from each sector, as described above, and 
then we implement the renormalization procedure of the 
temperature independent contributions. 

In the weak field limit, the contribution from the Higgs 
sector is given by 



V, 
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KUk2 
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(69) 



where the contributions from all scalars have been ac- 
counted for. 

The first term in Eq. ([M]) represents the lowest order 
contribution to the EP at finite temperature and zero 
external magnetic field, usually referred to as the boson 
ideal gas contribution |39| . The thermal piece of this 
contribution is given by (30[ 
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£f)+Oimf)y (70) 



327r2 V47rT 

Notice that there are potentially dangerous terms in 
Eq. (|70|) that can become imaginary for negative val- 
ues of mi. However as we will show, these terms cancel 
when including the Higgs contribution from the vacuum 
renormalization as well as the ring diagrams. 

The second, S-dependent term in Eq. (|69p vanishes 
identically [33 • Therefore, to one- loop order, the ther- 
mal contribution to the EP in the weak field case from 

I 



the Higgs sector is independent of eB and is given by 
Eq. dZOD. 

In the weak field limit, the contribution from the 
fermion sector is given by 



V, 



(1) 



N, 



2{eB) 






oJn + kl + mt'^ 



Y 



(71) 



where Nc = 3 is the number of colors. We emphasize 
that the only fermion mass we keep in the analysis is the 
top mass mt- 

The first term in Eq. (|7ip represents the fermion ideal 
gas contribution [39|, whose thermal part is explicitly 
given by 
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(72) 



The second term in Eq. (|7T|) is subdominant, after tak- 
ing care of renormalization as shown in Ref. 12J] . There- 
fore, to one-loop order, the thermal contribution to the 
EP in the weak field limit from the fermion sector is only 
given by Eq. ((7^ . 

In the gauge boson and ghost sectors, first note that 
in the broken phase the charged fields do couple to the 
external magnetic field. Therefore, their contribution to 
the EP at one-loop contains both magnetic field inde- 
pendent and dependent terms. For the former, the ghost 
field contributions cancel the spurious degrees of free- 
dom arising in the gauge sector when working in an ar- 
bitrary covariant gauge. In the latter, the charged ghost 
fields contribute as regular scalar fields with the oppo- 
site sign and vanishes identically [33[ [see discussion af- 
ter Eq. (|70| ]. Altogether, the magnetic field dependent 
contribution from these sectors comes only from the mag- 
netic field dependent pieces of the W boson propagator. 

Therefore, the thermal part of the contribution to the 
EP from the gauge boson sector is given explicitly by 
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where we have maintained the leading magnetic field de- 
pendent term and we defined 
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(74) 



Note that the factors in front of each magnetic field inde- 
pendent contribution in Eq. (|73p correspond to the two 
W's, the Z and the photon polarizations [23] ■ Also, no- 
tice that in the terms that come within the contributions 
proportional to the magnetic field squared, we have re- 
placed the Ty-boson mass by its thermal one within the 
argument of the logarithmic function and in the term 
where that mass appears in the denominator. This is be- 
cause, as previously discussed, although the analysis is 
strictly speaking only valid in the broken phase, since we 
consider the effective potential as a function of w > 0, to 
avoid divergences near w ~ we replace the gauge boson 
mass by the thermal one. 

To implement the renormalization procedure, notice 
that the T = one-loop EP from Eq. (p5)) involves the 
integral 



V(m) 



^2^ 



Vk^ 



(75) 



where m stands for the mass of any one of the contribut- 
ing species. This integral diverges and needs to be regu- 
larized. We do so by means of introducing the ultraviolet 
cutoff A and thus 
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(76) 



Therefore, the T = part of the one-loop potential can 
be written as 

3 

y{i)T=o ^ y (^^. A) + ^ V{m,; A) + 8V(mM/; A) 

3 

+ 4V(m2; A) - 12V(mt; A) - ^ V(7«^.; A), 

i=l 

(77) 

where the factors in front of each term account for the 
degrees of freedom of the corresponding species and the 
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FIG. 3. Renormalized one-loop effective potential for ^ = 1 
and different values of T, subtracting v independent terms 
and dividing by the fourth power of the t)-dependent W mass 
evaluated &\, v — v^ , the value of the vacuum expectation 
value that minimizes de effective potential after the phase 
transition is completed for a temperature T — 110 GeV. In 
this gauge the phase transition occurs at Tc — 128.9 GeV and 
is second order. 



masses depend on v. Notice that since we are working 
in arbitrary covariant gauge the unphysical degrees of 
freedom for gauge bosons are canceled by the contribu- 
tions of the ghost fields, which leaves Eq. ([77| with only 
the physical degrees of freedom. We emphasize that the 
only fermion we consider is the top quark. Since the 
theory is renormalizablc, it should be possible to absorb 
the A-dependent terms with the introduction of suitable 
counterternrs that maintain the form of the three-level 
potential. Therefore, the general structure of the T = 
effective potential up to one-loop order, after renormal- 
ization, should be 
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(78) 



where the coefficients a(A) and 6(A) are introduced to 
cancel the A-dependent terms in l/(i)^=" and the coeffi- 
cients 5c^ and 5\ take care of possible finite corrections of 
the z;2 and v^ terms, respectively. As the renormalization 
condition, we require that the minimum of Vren (v) re- 
mains at its classical value, namely 



UVren 
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(79) 



Introducing suitable constant terms to make the argu- 
ments of the logarithms be dimensionless and disregard- 
ing additive constants and terms explicitly proportional 
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to A^, we get 
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Therefore, after renormalization, the 1-loop EP at finite 
temperature is given by 
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where we have kept only the leading magnetic filed de- 
pendent term, with V2 as defined in Eq. (j74p and the 
masses without an argument are the u-dependent masses. 
As pointed out, the dangerous dependence of the scalar 
masses in the arguments of the logarithmic functions have 
disappeared. However, to cancel the cubic terms in the 
scalar masses one needs to go to the next order, namely, 
to the ring diagrams. Let us for the time being ignore 
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FIG. 4. The renomialized one-loop effective potential normal- 
ized by the fourth power of the v dependent W-mass evalu- 
ated at « = Uq', the value of the vacuum expectation value 
that minimizes de effective potential after the phase transition 
is completed for a temperature T = 110 GeV and C = li ^-s a 
function of the vacuum expectation value v and for different 
values of the gauge parameter ^. At this temperature, the 
position of the minimum wo shows an somewhat significant 
dependence on the gauge parameter for values ^ around 1. 



the terms proportional to cubic scalar masses and explore 
the properties of the renormalized effective potential at 

1-loop. Figure[3]shows Vrcn (w) for ^ = 1 and different val- 
ues of T, subtracting v independent terms and dividing 
by the fourth power of the f-dependent W mass evalu- 
ated at w = W(^, the value of the vacuum expectation value 
that minimizes de effective potential after the phase tran- 
sition is completed for a temperature T = 110 GeV. We 
can see that the phase transition is second order and in 
this gauge it happens for T^. = 128.9 GcV. To illustrate 
the gauge dependence of the phase transition parame- 
ters, Fig.|3]shows Vrcn (w) for a low temperature T = 110 
GcV, for which the phase transition has been completed, 
for different values of the gauge parameter ^. Notice the 
dependence on the gauge parameter of the broken phase 
minimum at this constant temperature, whose numeri- 
cal values are shown on the second column of table I. 
The dependence on the gauge parameter of the critical 
temperature is milder, as shown on the first column of 
tabic I. 



B. Ring diagrams 

It is well known that the next order correction to the 
EP comes from the so called ring diagrams. These are 
schematically depicted in Fig. [S] Their renormalized con- 
tribution to the EP is more clearly found by explicitly 
separating the two-loop contribution [39[ coming from 



14 



C Tc (GeV) vo (GeV) 

3.5 120.8 102.3 

1 128.9 130.9 

0.1 132.4 139.0 

TABLE I. Critical temperature Tc (second column) and posi- 
tion of the minimum vo for the constant temperature T = 110 
GeV (third column) of the renormalized one- loop effective po- 
tential, after the phase transition has been completed. The 
chosen values of ^ are close to 1, Notice that the critical tem- 
perature is only mildly dependent on 5, however the position 
of the minimum for a temperature where the phase transition 
has been completed shows a larger ^ dependence. 



the scalar sector and can be written as 
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2 ^ j (27r)3 
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(82) 



where H^, 11"-^ and n*^ are the Higgs, neutral and 
charged scalar self-energies in the presence of the mag- 
netic field, D^ , D^ and D'~^ their corresponding propa- 

(2) 

gators and V"^ the contribution to the two-loop effective 
potential coming exclusively from the scalar sector. The 
factor 2 accounts for the two charged scalar degrees of 
freedom. As written, Eq. ((82|) deserves some comments: 
First, the dominant contribution comes from the Matsub- 
ara frequency with n = 0. Second, ultraviolet divergences 
are canceled cxplicity. Third, the full two-loop contribu- 
tion involving the scalar sector contains diagrams with 
particles other than scalars. For the purposes of the 
present analysis, we consider only the subclass of dia- 
grams that contain only scalars and thus we get 



K 



(2) 



A 7^4 

24 



IGtt 
A 



1927rmJ 






TO3 + 2?71l) T 



(83) 



where we keep only the leading contributions, within the 
hierarchy of energy scales considered. The potentially 
dangerous terms with odd powers of the scalar masses 
in Eq. (|83p exactly cancel similar terms coming form the 
integral in Eq. (j82p . under the approximation that the 
self-energies involve only the scalar contribution. The 
full proof of this cancellation in the absence of magnetic 
fields has been treated in detail for the standard model in 
Ref. [2^. In the presence of an external magnetic field, 
it has been shown in the linear sigma model in Ref. [46| . 
Within our approximation, this result shows that con- 
sidering only scalar contributions, this cancellation also 



FIG. 5. Schematic representation of the resummation of the 
ring diagrams. 



happens in the standard model with a magnetic field. A 
full proof is currently under way [47 1. 

The dominant contribution in Eq. (|82p comes from the 
mode n = 0. The explicit expression for Eq. ()82p is [33] 



^ .(ring) _ '^ ^4 _ T r , 

^" "24 12^ U 
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As anticipated, the cubic mass terms in Eq. ((M)) will 
exactly cancel the similar terms appearing in Eq. (|8ip 
after adding both equations. 

Next, we turn to the ring contribution from the gauge 
boson and Faddeev-Popov ghost sectors. The compu- 
tation is simplified by recalling that the ghost degrees 
of freedom cancel the leading temperature contributions 
from spurious degrees of freedom arising form the covari- 
ant treatment of the gauge fields. The contribution that 
depends on the magnetic field, being subleading, is not 
canceled, as was the case for the one-loop contribution. 
Thus we write 



K 



(ring) 



{v) + V, 



(ring) 
FP 



T 



EE 



(pk 

1^ 



X Tr 






pAgb(0)G^a^«,k)] 



N 



y 



{ring) 



(85) 



charged FP 

where the sum over gb runs over the four gauge boson 

physical degrees of freedom and V^hartLd fp represents 
the subleading, magnetic field dependent term, arising 
from the ring contribution of the charged ghost fields. 
Notice that this time, in contrast to the scalar case, we 
do not add explicitly the two-loop contribution and con- 
sequently the sum over the index N starts from N = 2. 
We first compute the contribution from the neutral 
gauge bosons to Eq. (|85p . In the approximation where 
the particle has a small momentum (the infrared limit), 
the Euclidean version of a neutral, gauge boson propaga- 
tor, for an arbitrary value of the gauge parameter ^ can 
be written as 
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where 
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By using Eqs. ([65|) and ([87| . it is easy to see that the 
product n^A(0)G''*"'(a;„,k) becomes 
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from where, considering the n = term and taking the 
trace, we can write Eq. ([55]) as 
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Notice that the plasma screening effects in Eq. ([89]) are 
naturaUy accounted for by the fact that the thermal mod- 
ification of the gauge boson mass mgb {rhz and rhj) ap- 
pear. Also, the nice cancellations that took place in the 
scalar sector between odd powers of the masses will not 
happen in this case when adding these terms to Ec^. ()8ip . 



This is because in the infrared limit there is no contribu- 
tion from the gauge boson transverse degrees of freedom, 
thus there is no match of the coefficients to produce the 
cancellation. This is a feature of the infrared limit we 
are considering. However, since the square of the gauge 
boson masses are never negative, these terms do not pose 
a problem. 

Next, we turn to the contribution of the charged 
gauge fields. In the weak field limit, the product 
n^A(0)G'^''(a;„,k) can be written as 
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Using Eq. (j90|) . and carrying out an expansion of the 
argument of the logarithm, we can explicitly write 
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Using this result into Eq. (|85|) . considering the n = 
term and taking the trace in Eq. (|85p we get 
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r 



In this approximation, the plasma screening effects in 
Eq. (|92p emerge naturally as a thermal modification of 
the gauge boson mass, rh-w defined in Eqs. (p^ and (|55| . 
Finally, the magnetic field dependent contribution 
from the charged ghost fields to the ring potential is eas- 
ily computed recalling that this contribution is equivalent 
to the one coming from two scalars (one for each W field) 
but with opposite sign. Thus we have 
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where the diagram depicting the Faddeev-Popov ghosts 
is shown in Fig. [5] and, as shown in Appendix C, the 



leading contribution to the W^-ghost self energy is given 
by 
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and the thermal iy(Z)-ghost mass is given by 
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S : Higgs and Goldstone bosons 
G : 7, Z and W gauge bosons 
Gh : ghosts 



FIG. 6. Feynman diagrams for the Faddeev-Popov ghosts 
self-energies. Notice that since we work in the infrared limit, 
the diagrams of type (a), do not contribute to the ghosts self- 
energies. 



mw and mz arc the w-dependent W and Z masses, re- 
spectively and "fhi {i = 1 ... 4) are defined in Eqs. (|3]) 
and (|17|) . Notice that in Eq. ([M)) . we have replaced rrii 
by rhi since, although the analysis is valid near the bro- 
ken phase minimuni, we consider the effective potential 
as a function of w > 0, and for small values of v, the 
square of the scalar masses can become negative. 



C. Effective potential up to ring order 



The final expression for the effective potential 
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is obtained by adding up the results in Eqs. ([55)1 . ([Tn]). 
dZl]), dZSl), dill), dMl) and (HSl). 

In order for the terms involving the square of the scalar 
bosons' thermal mass to be real, the temperature must 
be such that 
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which defines a lower bound for the temperature. A more 
restrictive bound is obtained by requiring that eB < mj^ 
for the weak field expansion to work. This condition 



translates into the bound 
T>T2 = 



eB + 16c2 



3g2 -f g'2 + 8A + 4/2 



(98) 



The relevant factor that enhances the order of transition. 



present both in V^^™ 



and V„ 



(ring) 



which 



.g, , is {eBY/mf 
can be traced back to the boson self-energy diagrams 
involving a tadpole of charged scalars in the presence of 
the external field. 



VI. SYMMETRY RESTORATION 

In order to quantitatively check the effect of the mag- 
netic field during the EWPT, we proceed to plot Voff as 
a function of the vacuum expectation value v. For the 
analysis we use g' = 0.344 and g = 0.637, mz = 91 GeV, 
■mw = 80 GeV, / = 1, A = 0.11 which corresponds to 
the current bound on the Higgs mass. 

Figure [7] shows the effective potential in the absence 
of magnetic field divided by [r7T.vi/(wf )]'' = (10.5 GeV)''^, 
where vj is the value where the broken phase minimum 
appears at the critical temperature, which in this case 
happens to be T = 139.758 GeV. The value of the gauge 
parameter is ^ = 0.1. The phase transition is weakly first 
order. Note that the results are consistent with the ones 
obtained in Rcf. [2J] which are computed using symmet- 
ric phase degrees of freedom. 

Figure [5] shows the effective potential divided by 
[?TivF('i^i")]^ = (10.5 GeV)^ for the same temperatures 
as in Fig. [7] and a fixed value of the magnetic field 
parametrized as B = 6 x (100 GeV)2, with b = 0.01. 
An insert in this plot shows the difference (AVg^) with 
respect to the effective potential shown in Fig.[7]m the ab- 
sence of magnetic field over a small region in the v range 
where the second minimum would be developing and is 
now delayed by the presence of magnetic field. This same 
effect can be observed if we keep the temperature fixed 
and increase the value of the magnetic field. Given that 
we use a weak magnetic field, the effect is small. In order 
to appreciate such effect, in Fig. [5] we show the difference 
between the effective potentials in the presence and in 
the absence of magnetic field. Starting from zero mag- 
netic field, for which the phase transition happens at the 
critical temperature T ~ 139.758 GeV, the phase transi- 
tion is delayed by increasing the values of the magnetic 
field, also parametrized as B = bx (100 GeV)2, while the 
temperature is maintained fixed. In both Figs. [S] and [HI 
the value of the gauge parameter is ^ = 0.1. 

Figure [To] shows a small region around the second min- 
ima that develops in the effective potential divided by 
the fourth power of the u-dependent W mass evaluated 
at vf , for three different values of the magnetic field and 
computed at their corresponding critical temperatures, 
keeping ^ = 0.1. Notice that for increasing values of the 
magnetic field, the phase transition starts at a lower crit- 
ical temperature Tc and the broken phase minimum vq is 
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FIG. 7. Effective potential divided by the fourth power of the 
w-dependent W mass evaluated at Vi (the value where the 
broken phase minimum appears at the critical temperature in 
the absence of magnetic field). The gauge parameter ^ = 0.1 
and the value of the magnetic field is set to zero. The phase 
transition is weakly first order. Note that the results are 
consistent with the ones obtained by [24| which are computed 
using symmetric phase degrees of freedom 
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FIG. 9. Difference between the effective potentials in the 
presence and in the absence of magnetic field divided by the 
fourth power of the u-dependent W mass evaluated at vi (the 
value where the broken phase minimum appears at the crit- 
ical temperature T = T^^" = 139.758 GeV in the absence 
of magnetic field) for different values of the magnetic field 
parametrized as B = b x (100 GeV)^. The gauge parameter 
is 5 = 0.1. Even though the magnetic field is weak, we can ap- 
preciate that the phase transition is delayed when increasing 
the values of the magnetic field. 



also shifted to higher values in such a way that the ratio 
vq/Tc increases at the phase transition. This increase, 
although modest, is a desired feature that can eventually 
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FIG. 8. Effective potential divided by the fourth power of the 
u-dependent W mass evaluated at v^ (the value where the 
broken phase minimum appears at the critical temperature in 
the absence of magnetic field) for the same temperatures as in 
Fig. [7] and a fixed value of the magnetic field parametrized as 
B = b X (100 GeV)^, with b = 0.01. The gauge parameter is 
5 = 0.1. The insert shows the difference (AVg^) with respect 
to the effective potential shown in Fig. \7\ in the absence of 
magnetic field over a small region in the v range where the 
second minimum would be developing and is now delayed by 
the presence of magnetic field. 



help a possible baryon asymmetry to not be washed out 
after the completion of the phase transition Q . 

To explore the gauge parameter dependence of the ef- 
fective potential's relevant parameters, Figs.fTTj-fT^show 
the difference in the behavior with and without magnetic 
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FIG. 10. A small region around the second minima that de- 
velops in the effective potential divided by the fourth power 
of the ^-dependent W mass evaluated at vf (the value where 
the broken phase minimum appears at the critical tempera- 
ture in the absence of magnetic field) for three different values 
of the magnetic field at their corresponding critical tempera- 
tures and for ^ = 0.1. Notice how for increasing values of the 
magnetic field, the phase transition starts at a lower temper- 
ature and the broken phase minimum is also shifted to higher 
values in such a way that the ratio vq/Tc increases at the 
phase transition. 
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FIG. 11. Difference of gauge parameter dependence of uq with 
and witliout magnetic field. Tire position of the broken phase 
minimum at the critical temperature, for two values of the 
magnetic field parametrized as _B = 6 x (100 GeV)^ remains 
stable throughout a large range of ^, even in the presence of 
magnetic field. 



FIG. 13. Difference of gauge parameter dependence of the 
ratio vo/Tc with and without magnetic field, for two values 
of the magnetic field parametrized as B = fe x (100 GeV)^. 
Notice how the value of vo/Tc is stable for a variation of ^ up 
to 1.5. The presence of the magnetic field does not introduce 
a strong gauge dependence. 



field oivQ, Tc and the ratio vq/Tc for values of ^ around 1. 
We use a fixed value of the magnetic field, parametrized 
asB = bx (100 GeV)2, with 6 = 0, 0.005, 0.01. We can 
see that wg, T^. and Vq/T^ remain stable for a variation 
of ^ up to 1.5. For these small values of the magnetic 
field strcnght, consistent with the assumed hierarchy of 
energy scales, the magnetic field does not introduce a 
strong gauge dependence. 
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FIG. 12. Difference of gauge parameter dependence of Tc 
with and without magnetic field. The position of the critical 
temperature for the development of the broken phase min- 
imum, for two values of the magnetic field parametrized as 
B — b X (100 GeV)^ remains stable throughout a large range 
ofC. 



VII. DISCUSSION AND CONCLUSIONS 

In conclusion, we have studied the symmetry restora- 
tion problem in the MSM at finite temperature in the 
presence of an external magnetic field. The study has 
been carried out by analyzing the finite temperature ef- 
fective potential. We have shown that in the presence 
of the magnetic field, the EWPT has become a stronger 
first order phase transition. Our treatment has been im- 
plemented for the case of weak magnetic fields for the 
hierarchy of scales eB ^ m^ <^ T^ where m is taken as 
a generic mass involved in the calculation. 

We have explicitly worked with the degrees of free- 
dom in the broken symmetry phase where the external 
magnetic field belongs to the U{l)em group and in a co- 
variant gauge, with an arbitrary value of a single gauge 
fixing parameter. The calculation is carried out up to 
the contribution of ring diagrams. To include the effects 
of the magnetic field, we have made use of the Schwinger 
proper-time method to describe the particle's propaga- 
tors. In this way, the contribution from all Landau levels 
has been accounted for. We have carried out a system- 
atic expansion up to order (eB)'^ and have computed, as 
an intermediate step, the magnetic field and gauge fixing 
parameter dependent self-energies in the MSM J48l |. 

The presence of the external magnetic field gives rise 
to terms in the effective potential proportional to 1/m?, 
where rhf = mf + IIi, coming from tadpole diagrams 
in the boson self-energies where the loop particle is a 
charged scalar and m,, are their masses. These terms are 
the relevant ones for the strengthening of the order of the 
phase transition. The results are in qualitative and quan- 
titative agreement with the ones previously found [2J] us- 
ing symmetry restored degrees of freedom. The ratio of 
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the broken phase vacuum expectation value to the criti- 
cal temperature at the phase transition increases with in- 
creasing values of the magnetic field. This feature works 
in favor of the suppression of the sphaleron induced tran- 
sitions in the broken phase after the EWPT is completed 
and thus against a washing out of a possibly produced 
baryon number. The increase in this ratio is modest, 
though it is important to notice that wc have used very 
restrictive values of the field strength, consistent with the 
working hierarchy of energy scales used in this work. The 



leading temperature terms of the effective potential are 
gauge parameter independent, although the subleading 
ones, i.e., the ones coming with the magnetic field turn 
out to be gauge parameter dependent. The dependence 
on this parameter of phase transition observables such 
as the position of the minimum in the symmetry broken 
phase and the critical temperature for the transition is 
not negligible for values of ^ near ^ ~ 1. This dependence 
signals that the truncation of the magnetic field power se- 
ries might not be a gauge invariant procedure. This point 
deserves a closer look and is being investigated |47| . 
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APPENDIX A: SCALAR SELF-ENERGIES IN ARBITRARY GAUGE 

In this appendix, we list the results for the self-energy diagrams depicted in Fig. [T] (for the neutral Goldstone boson 
there is no contribution of the type shown in diagram (c) of Fig. [T]). 



A. Higgs boson {H) self-energy in arbitrary gauge 
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Neutral Goldstone boson (G") self-energy in arbitrary gauge 
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Charged Goldstone boson (G"^) self-energy in arbitrary gauge 
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APPENDIX B: GAUGE BOSON SELF-ENERGIES IN ARBITRARY GAUGE 

In this appendix, we list the resuhs for the self-energy diagrams depicted in Fig. [5] 



Photon (7) self-energy in arbitrary gauge 
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Neutral gauge boson (Z) self-energy in arbitrary gauge 
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Charged gauge boson (W) self-energy in arbitrary gauge 
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APPENDIX C: GHOSTS SELF-ENERGIES IN ARBITRARY GAUGE 



In this appendix, we list the resuhs for the self-energy diagrams depicted in Fig. [6l 



7 boson associated ghost (r;^) self-energy in arbitrary gauge 
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Z boson associated ghost (772) self-energy in arbitrary gauge 
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W boson associated ghost (riw) self-energy in arbitrary gauge 
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